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CN ; Abstract 
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• One-dimensional system of Brownian motions called Dyson's model is the 

Z/y . particle system with long-range repulsive forces acting between any pair of 

I particles, where the strength of force is P/2 times the inverse of particle 

' distance. When /? = 2, it is realized as the Brownian motions in one di- 

,— 1< mension conditioned never to collide with each other. For any initial con- 

I figuration, it is proved that Dyson's model with (3 = 2 and N particles. 

Oh; X{t) = {Xi{t),. . . ,XN{t)),t G [0,00), 2 < AT < 00, is determinantal in 

■ the sense that any multitime correlation function is given by a determinant 

, with a continuous kernel. The Airy function Ai(z) is an entire function 

^ I with zeros all located on the negative part of the real axis M. We con- 

sider Dyson's model with (3 = 2 starting from the first N zeros of Ai(z), 
t:;^- . 0>ai>--->aAr,A>2. In order to properly control the effect of such ini- 

\ tial confinement of particles in the negative region of M, we put the drift term 

^5 ■ to each Brownian motion, which increases in time as a parabolic function : 

\0 '. Yj{t) = Xj{t)+ty4+{di+J2f=i{l/ae)}t, l<j<N, where di = Ai'(0)/Ai(0). 

; We show that, as the TV ^ 00 limit of Y{t) = (Yiit), . . . ,YNit)),t e [0,oo), 

' we obtain an infinite particle system, which is the relaxation process from the 

. configuration, in which every zero of Ai(2;) on the negative M is occupied by 

Q [ one particle, to the stationary state /iAi- The stationary state /iAi is the deter- 

' minantal point process with the Airy kernel, which is spatially inhomogeneous 

on M and in which the Tracy- Widom distribution describes the rightmost par- 
ticle position. 
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1 Introduction 



1.1 Dyson's model: One-Dimensional Brownian Particle Sys- 
tem Interacting through Pair Force \ jx 

To understand the time-evolution of distributions of interacting particle systems on 
a large space-time scale [thermodynamic and hydrodynamic limits) is one of the main 
topics of statistical physics. If the interactions among particles are short ranged, 
the standard theory is useful. If they are long ranged, however, general theory has 
not yet been established and thus detailed study of model systems is required [26] . 

In the present paper, we consider Brownian particles in one dimension with long- 
ranged repulsive forces acting between any pair of particles, where the strength of 
force is exactly equal to 1/x when the particle distance is x. If the number of particles 
is finite N < oo, the system is described by E(t) = ^Xj{t), 2 < N < oo, where 
X{t) = (Xi(t), . . . , XAr(t)) satisfies the following system of stochastic differential 
equations (SDEs); 

dX^it) = dB,{t) + J2 x(t)-X (t) "^^' 1<^'<^' te[0,oo) (1.1) 

with independent one-dimensional standard Brownian motions Bj{t),l < j < N. 
The SDEs obtained by replacing the 1/x force in (11. ip by /3/{2x) with a parameter 
P > were introduced by Dyson |1] to understand the statistics of eigenvalues of 
hermitian random matrices as particle distributions of interacting Brownian mo- 
tions in M. Corresponding to the special values (3 = 1,2 and 4, hermitian random 
matrices are in the three statistical ensembles with different symmetries, called the 
Gaussian orthogonal ensemble (GOE), the Gaussian unitary ensemble (GUE), and 
the Gaussian symplectic ensemble (GSE), respectively [3]. In particular for /? = 2, 
that is the case of (II. ip . if the eigenvalue distribution of N x N hermitian random 
matrices in the GUE with variance cr^ is denoted by /i^^f , we can show 

hm = f^sU-): (1-2) 



iV- 



where /ism denotes the determinantal (Fermion) point process pll [23 j with the so- 
called sine kernel 

KUx) = —f dke^"'^ =^-^^^^, xeR. (1.3) 
271 y|fc|<^ nx 

That is, fisin is a spatially homogeneous particle distribution, in which the particle 
density is given by psin = linix^o -^sin(a;) = 1 and any A^i-point correlation function 
Psm{xN-^),XN-^ = {xi, . . . ,xni) G M.'^'^ , Ni > 2, is given by a determinant of an 
Ni X Ni real symmetric matrix; 



PsinixNi) = det 

l<j,k<Ni 
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Based on this fact known in the random matrix theory Spohn |25j studied 
the equihbrium dynamics obtained in the infinite-particle hmit iV — >^ cx) of Dyson's 
model fll.ll) . Since the Xjx force is not summable, in the infinite-particle limit 
N ^ oo the sum in f ll.ll) should be regarded as an improper sum, in the sense that 
for Xj{t) G [—L,L] the summation is restricted to fc's such that Xkit) G [—L,L] 
and then the limit L — oo is taken. It is expected that the dynamics with an 
infinite number of particles can exist only for initial configurations having the same 
asymptotic density to the right and left [IS]- 

The problem, which we address in the present paper, is how we can control 
Dyson's model with an infinite number of particles starting from asymmetric initial 
configurations. The motivation is again coming from the random matrix theory as 
follows. Consider the Airy function [H [30] 

Ai(z) = — / dke^^'^+'^'/'^l (1.4) 

It is a solution of Airy's equation f"{z) — zf{z) = with the asymptotics on the 
real axis M: 

Ai(x)^^-^exp (-^x^/^ 

Ai(— x) ~ — ^^^^ cos ^-x'^/^ — — ^ in x ^ +oo. (1.5) 

In the GUE random matrix theory, the following scaling limit has been extensively 
studied: 

lim /i^u%(2iv2/3 + .)=^^^(.), (1.6) 

N—^00 ' 

where /^Ai is the determinantal point process such that the correlation kernel is given 
by PI2H], 

du Ai{u + x)Ai(M + y) 

^ Ai(x)Ai-(|/)-Ai-(x)Ai(y) 

^y ' ^^^^^ (1.7) 

(Ai'(x))2 - x(Ai(x))2, x = y eR. 

It is another infinite-particle limit different from (11.21) and is called the soft-edge scal- 
ing limit, since x^/2t ^ (2A^^/^)^/(2A^^/^) = 2A^ marks the right edge of semicircle- 
shaped profile of the GUE eigenvalue distribution (see, for example, [12] )• The 
particle distribution fiAi with the Airy kernel (11.71) is highly asymmetric: As a mat- 
ter of fact, the particle density pAi{,x) = KAi{x\x) decays rapidly to zero as x — >• oo, 
but it diverges 

PAi{x) ~ — (— x)"*^^^ ^ oo as X ^ — oo. (1.8) 

TT 
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Let R be the position of the rightmost particle on M in yUAi- Then its distribution is 
given by the celebrated Tracy- Widom distribution [2H] 



/iAi(-R < x) = exp 



{y ~ x){q{y)fdy 



where q{x) is the unique solution of the Painleve II equation 

// I o 3 

q = xq + 2q 

satisfying the boundary condition q{x) ~ Ai(x) in x — > oo. Prahofer and Spohn 
and Johansson [TT] studied the equilibrium fluctuation of this rightmost particle and 
called it the Airy process. Tracy and Widom derived a system of partial differential 
equations, which govern the Airy process [29] . See also [21 [3] . How can we realize 
fiAi as the equilibrium state of Brownian infinite-particle system interacting through 
pair force 1/x ? The initial configurations should be asymmetric, but what kinds 
of conditions should be satisfied by them ? How should we modify the SDEs of 
original Dyson's model fll.ll) . when we provide finite-particle approximations for 
such asymmetric infinite particle systems ? 

In the present paper, as an explicit answer to the above questions, we will present 
a relaxation process with an infinite number of particles converging to the stationary 
state /iAi in t ^ oo. Its initial configuration is given by 

oo 

U-) = Y.^ai-) = Y.^a,i-), (1.9) 

a€A j=l 

in which every zero of the Airy function fll.4p is occupied by one particle. This 
special choice of the initial configuration is due to the fact that the zeros of the Airy 
function are located only on the negative part of the real axis M, 

A = Ar\0) = joj, J G N : Ai(aj) = 0, > ai > > ■ ■ ■ |, (1.10) 

with the values [I] Oi = —2.33 . . . , a2 = —4.08 . . . , 03 = —5.52 . . . , 04 = —6.78 . . . , 
and that they admit the asymptotics [H 130] 

— in 7^00. (1.11) 



«i - - I Y ) J ' m J 

Then the average density of zeros of the Airy function around x, denoted by PAi-i(o) (^) ) 
behaves as ^ 

PAi-i(o)(^) — ^00 as X ^ —00, 

which coincides with (11.81) . The approximation of our process with a finite number 
of particles < 00 is given by S^(t) = XljLi ^Vjit) with 

r,(t) = x,(t) + ^ + D^^t, i<j<N, te[o,oo), (1.12) 
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associated with the solution Xit) = {Xi{t), . . . , Xjyit)) of Dyson's model (11.11) . 
where 

^ 1 

Here di = Ai'(0)/Ai(0) and An = |o > ai > • ■ ■ > ajy^ C ^ is the sequence of the 

first zeros of the Airy function. In other words, Y{t) = {Yi{t), ¥2(1), ... , YN{t)) 
satisfies the following SDEs ; 



l<fc<iV 

ky^j 



l<j<N, te[0,oo), (1.14) 



where Bj{t)'s are independent one- dimensional standard Brownian motions. 

For 1^(0) = X E M^, set ^^(■) = J2f=i^xj{-) and consider the process S^(t) 
starting from the configuration . We consider a set of initial configurations 
such that they are in general different from the A^-particle approximation of (11.91) . 

AT 

^a{-)= E ^^(■) = T.^^A-)^ (1-15) 

aeAM j=i 

but the particle density p{x) of limAr_^oo^^ will show the same asymptotic in x — > 
—00 as (II. 8p . Because of the strong repulsive forces acting between particle pairs in 
(II. ip . such confinement of particles in the negative region of M at the initial time 
causes strong positive drifts of Brownian particles. The coefficient (I1.13P of the drift 
term -D^^t added in (I1.12p . however, negatively diverges 



1/3 



/12\ 

(;^) A^^/^^-00 as A^^oo. (1.16) 



We will determine a class of asymmetric initial configurations denoted by Xff, which 
includes as a typical one, such that the effect on dynamics of asymmetry in 
configuration will be compensated by the additional drift term D_A^t in the infinite- 
particle limit N —>■ 00 and the dynamics with an infinite number of particles exists. 
Note that we should take N 00 limit for finite t < cxd in our process (11.120 to 
discuss non-equilibrium dynamics with an infinite number of particles. In the class 
X^, when the initial configuration is specially set to be (II. 9p . we can prove that 
the dynamics shows a relaxation in the long-term limit t — 00 to the equilibrium 
dynamics in //Ai- In the proof we use the special property of the systems (II. ip 
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and fll.l4p such that the processes have space-time determinantal correlations. This 
feature comes from the fact that if and only if the strength of pair force is exactly 
equal to 1/a; when the particle distance is x, i.e., iS jS = 2, Dyson's model is realized 
as the Brownian motions conditioned never to collide with each other [TOt [T2]. 

In order to explain the importance of the notion of entire functions for the present 
problem, we rewrite the results reported in our previous paper [13j for Dyson's model 
with symmetric initial configurations below. Then the changes which we have to do 
for the systems with asymmetric initial configurations are shown. There the origin 
of the quadratic term in (11.121) will be clarified. 

1.2 Processes with Space-time Determinantal Correlations 
and Entire Functions 



In an earlier paper [IB], we studied a class of a non-equilibrium dynamics of Dyson's 
model with {3 = 2 and an infinite number of particles. As an example in the class, 
we reported a relaxation process, denoted here by (S(t),Psin), which starts from a 
configuration 

U-) = Y.^a{-), (1.17) 

in which every point of Z is occupied by one particle, and converges to the sta- 
tionary state /isin. This process (S(t),Psin) is determinantal, in the sense that 
there is a function Ksin(s, x; t, ?/) called the correlation kernel such that it is con- 
tinuous with respect to {x,y) G for any fixed (s,t) G [0, oo)^, and that, for 
any integer M > 1, any sequence {N^)^^^ of positive integers, and any time se- 
quence < ti < ■ • • < tftf < oo, the {Ni, . . . , NM)-multitime correlation function 
Psm(ii, a^SiJ; . . . ; tM, x^^J), x^^J^ = {x^^\ xj^^) G M^™, 1 < m < M, is expressed 
by a determinant of a J2m=i x X]m=i asymmetric real matrix; 



L [til X]^_^] ■ ■ ■ ', tjyj , Xj^^^ j — det lKsin(im5 Xj ] tji, X 

l<m,n<M 



:i.i8) 



The finite dimensional distributions of the process (S(t),Psin) are determined by 
Ksin through (11.181) . It is expected that the correlation kernel Ksin is described by 
using the sine function as is the correlation kernel i^^gm of the stationary distribution 
/^sin given by (II. 3p . It is indeed true. Set 



and 



/(z) = sin(7r;2), z G C, (1.19) 

V,iJt,x) = ^= tGM\{0}, xGC. (1.20) 
y 27r|t| 

When t > 0, psm{t,y — x) is the heat kernel: the solution of the heat equation 
du{t,x)/dt = {l/2)d'^u{t,x) /dx"^ with \mit^Qu{t,x)dx = Sy{dx), and is expressed 
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using fll.lQp as 

pUt.y- ^) = \j^due--''^''/^[f{ux)f{uy) + f{ux + 1/2) f{uy + 1/2) 
For < s < t, by setting (11.201) . the Chapman- Kolmogorov equation 

/ dypsin{t - S,Z - y)psin{s,y - X) = Psin{t,Z - x) (1.21) 

Jr 

can be extended to 

dypsm{-t, z - y)psm{t - s,y-x) = Psin(-s, z-x). (1.22) 

Then ]Ksin(s, x; t, y) is given by 

, f dz 1 f(z) 

Kf{s,x;t,y) = > / ^=p/(0, a; s, x) -—-pf{t,y;0,z) 

naf-un^J^/^^^-^ z-a/(aj 



ae/-i(0)' 

-l{s > t)pf{t,y;s,x), s,t>0, x,y eR (1.23) 

with setting (I1.19P and pf{s,x]t,y) = Psm(t — s,y — x) with (ll.20p . where f~^{0) 
denotes the zero sei of the function /; /~^(0) = {z : f{z) = 0}, f'{a) = df{z)/dz\z=a, 
and l(c<j) is the indicator of a condition u; 1{uj) = 1 if is satisfied, and 1{uj) = 
otherwise. In this paper dz ■ means the integral on the imaginary axis in C 

from —\/—loo to a/— loo. The well-definedness of the correlation kernel and 
thus of the process (S(t),Psin) is guaranteed [I3] by the fact that the sine function 
(I1.19P is an entire function {i.e., analytic in the whole complex plane C), and the 
order of growth pj, which is generally defined for an entire function / by 

loglogM/(r) 

p/ = limsup for iW/(r) = max 1 7(2;) I, 

r^oo log r \z\=r 

is one. (The type defined by o"/ = limsup^^o^ log M/(r)/r''-^ is equal to n for (11.191) . 
That is, the sine function (I1.19P is an entire function of exponential type vr [H] ; 
Mf{r) ~ e'"^ as r — s> 00.) We can show that 

K^init,x;t,y)Ksinit,y;'t,x)dxdy ^ ^z{dx)l{x = y) as t ^ 0, 

since /~^(0) = sin~^(0)/7r = Z. It implies that the initial configuration (I1.17P of 
the relaxation process (S(i:),Psin) shall be regarded as the point-mass distribution 
on the zero set of the sine function (11.190 . Moreover, we showed in [13], by noting 



1 m 

-af'ia) 



Ksin{z - a), 
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if a G f-\0) =Z and z^a for ffLTOj) . that 

^sin{s + 9,x]t + e,y) ^ Ksin{t - s,y - x) as 6*^00 
with the so-called extended sine kernel, 



(1.24) 



du e''^"^*/^ cos(7rMx) 



if t > 
if t = 



POO 

- du e"'"'*/2 cos(7rMa;) if t < 
Ji 



[1.25) 



X e M. The equilibrium dynamics in yUgin, first studied by Spohn [25], has been shown 
to be determinantal with the correlation kernel fll.25p by Nagao and Forrester [18]. 
This process is realized in the long-term limit of the relaxation process Psin)- 
See also the Dirichlet form approach by Osada to the reversible process with respect 

to [201 EI]. 

Now we set 

f{z) = Ai{z), zeC. (1.26) 

The Airy function Ai{z), (11. 4p . is another entire function, whose order of growth is 
Pf = 3/2 with type aj = 2/3; max|^|=r \fiz)\ ~ exp[(2/3)r^/^] as r ^ 00. For t > 0, 
we consider 

PAiit,y\x)= [ due^'^^f{u + x)f{u + y), x,i/GM, (1.27) 

which is the solution of the differential equation; 
d 



dt 



1/92 

u(t,x) = - { -—r — X ] u(t,x) with \imu(t,x)dx = Sy(dx). 
2 V ox^ J t^o 



The integral / dzpAi{t, z\x) is given by 



g{t,x) = exp ( -y + — 



:i.28) 



We find, for s < t < 0, g{s,x)pAi{t — s,y\x)/g{t,y) is equal to the transition proba- 
bility density of 



Bit) + - 



;i.29) 



from X at time s to ?/ at time t, x,y G M, where B{t),t G [0, 00) is the one- 
dimensional standard Brownian motion. (See also [6J and references therein.) Then 
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for s,t E M., s ^ t, x,y E C, we set 



q{s, t,y-x) = psin { t - s, { y - — \ - \ X 



exp 



{y-xf , {t + s){y-x) {t-s){t + sy 
2{t - s) 



+ 



4 



32 



;i.3o) 



and as an extension of fll.27p we define 

PAi{t - s, y\x) = — -q{s, t,y-x) 



g{s,x) 



v/27r|t 



exp 



(y-xf {t-s){y + x) , (t-sf 
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2(t - s) 4 
Corresponding to fll.211) and (11.221) . we have the two sets of equahties 

dyq{s,t,z-y)q{0,s,y-x) = q{0,t,z-x) 



and 



/ 



dyq{t,0,z - y)q{s,t,y - x) = q{s,0,z-x) 

dypAiit - s, z\y)pAiis,y\x) = pAiit, z\x) 
dypAii-t, z\y)pAi{t - s,y\x) = PAi{-s, z\x) 



1.31) 



^1.32) 
:i.33) 

^1.34) 
:i.35) 



for < s < t. 

Let be the function given by (11.231) with setting (I1.26P and pf{s,x;t,y) = 
PAiit — s,y\x) with (ll.3ip . We will prove that Kai is well-defined as a correlation 
kernel and it determines finite dimensional distributions of an infinite particle system 
through a similar formula to (11.180 . We denote this system by (S_4(t), Pai)- The 
fact that pAi used in KAi is a transform (11.310 of the transition probability density 
q of (11.290 is the origin of the quadratic term t^/4 in (I1.12p . We can show that 



KAiit,x;t,y)KAi(t,y]t,x)dxdy ^ C,j[{dx)l{x = y) as t 

By using the integral formula for (I1.26P 

1 f(z) 1 

^ ' ' du Ai(u + z)A\{u + a) 



0. 



z-af\a) (Ai'(a))2 7o 

for a E A, z ^ a, and the fact that {Ai(x + a)/Ai'(a), a E A\ forms a complete or- 
thonormal basis for the space 1/^(0, oo) of square integrable functions on the interval 
(0, oo) [2Z], we will prove that 



KAi(s + 6*, x; t + 6^, y) — ^ KAi(t — s, y|x) as oo, 



;i.36) 
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where Kai is the so-called extended Airy kernel^ 

< poo 

/ du e""*/^Ai(M + x)Ai{u + y) if t > 
KAi{t,y\x)=l ^ (1.37) 
-/ due""*/2Ai(u + x)Ai(u + i/) ift<0, 

^ J — oo 

x,y E M. We denote by (S_4(t),PAi) the infinite particle system, which is deter- 
minantal with the correlation kernel Kai [SI [IHl [ISl 12] ■ The Airy kernel (11.71) of 
/iAi is given by KAi{y\x) = KAi(0,?/|x) and thus (Syi(t),PAi) is a reversible pro- 
cess with respect to /xai- The process (S^(t),PAi), which is determinantal with the 
correlation kernel KAi, is a non-equilibrium infinite particle system exhibiting the 
relaxation phenomenon from the initial configuration to the stationary state /lAi- 

\ be the distri- 
bution of the process H^(t) = ^jLi5yj(i) starting from a configuration . We 
denote by Wl the space of nonnegative integer-valued Radon measures on M, which 
is a Polish space with the vague topology: we say C,n converges to ^ vaguely, if 

lim / (p{x)^nidx) = / (p{x)^{dx) for any (p G Co(M), where Co(M) is the set of all 



continuous real-valued functions with compact supports. Any element of 971 can 
be represented as ^(■) = J2jeA^^j(') with an index set A and a sequence of points 
in M, = {xj)j^\ satisfying ,^(/) = \l\{xj : Xj E 1} < oo for any compact subset 

/ C M. For A C M, we write the restriction of ^ on A as {C, fl A){-) = ^Xj{')- 

We put Mo = 1^ e OJl : ^({x}) < 1 for any x G m|. We will prove that the finite 

particle process (S_4(t),P^ ) is determinantal for any initial configuration G DJIq 
and give the correlation kernel Kj^ (Proposition 12. 4p . For G OK with an infinite 
number of particles ^(M) = oo, when ]K^'~^'^' converges to a continuous function as 
L — > oo, the limit is written as K^. If p^'"'^''^] converges to a probability measure 
on 1X)T[°'°°-*, which is determinantal with the correlation kernel K^, weakly in the 
sense of finite dimensional distributions as L ^ cxd in the vague topology, we say 
that the process (S^(t),P^) is well defined with the correlation kernel K^. (The 
regularity of the sample paths of H_4(t) will be discussed in the forthcoming paper 
|14j.) We will give sufficient conditions for initial configurations ^ G DJIq so that 
the process (S_4(t),P^) is well defined (Theorem 12.51) . We denote by X-^ the set 
of configurations ^ satisfying the conditions and put = X"^ fl DJIq. It is clear 
that the configuration G Xq-. Then, if we consider the finite particle systems 
Syi(^) = Ylf=i^Yj{t), N > 2, with (I1.12P starting from the A^-particle approximation 

of (ll.lSp . we can prove (H_4(t), P_^) (S_4(t), PAi) as ^ oo in the sense of fi- 
nite dimensional distributions (Theorem l2.6l (i)). That is, (S_4(t),PAi) = (H_4(t),P^) 
with (II. 9p . Moreover, we will show (I1.36P and prove the relaxation phenomenon 
(S^(t + 0),PAi) ^ (S^(t),PAi) (Theorem ESKii)). 



10 



The paper is organized as follows. In Sect. 2 preliminaries and main results 
are given. Some remarks on extensions of the present results are also given there. 
In Sect. 3 the properties of the Airy function used in this paper are summarized. 
Section 4 is devoted to proofs of results. 



2 Preliminaries and Main Results 

For ^(•) = X^jeA (■) £ we introduce the following operations; 
(shift) for ueR, T„e(-) = 

(square) ^^'\-) = J2^4-). 

We use the convention such that 

Ylf{x)^expl U{dx)\ogf{x)\= H 

for ^ G 971 and a function / on M, where supp ^ = {x G M : ^({a;}) > 0}. For a 
multivariate symmetric function g we write g{{x)xe£.) for g{{xj)j^\). 



2.1 Determinantal processes 

As an 9Jt-valued process P^), we consider the system such that, for any integer 

M >1, fm e Co(M), 6^^ G M, 1 < m < M, < ti < • • • < tju < oo, the expectation 

M „ 

of exp < 9m / fm{x)'E{trm dx) \ Can be expanded by 

^ m=l ^ 



Xr 



as 



exp< y^drn / fm{x)'E{tm,dx) 

U=i ^» 

l[dx 



iVi>0 NM>Om=l 

M N, 

X 

m=l j=l 



(1) 



Yldx 



(M) 



(M)\ 



(2.1) 



Here p's are locally integrable functions, which are symmetric in the sense that 
p(. . . ; tm, a(a;M); . . . ) = p(. . . ; t^, x^; . . . ) with a(a;(-)) = (0:^7,) , . . . , x^^^lj for 
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any permutation a G Sn^i 1 ^ ^ such a system p(ti, x^-*^); . . . ; a?*-^^-*) 

is called the (iVi, . . . ^ N M)-fnultitime correlation function and Q^{x\ the generating 
function of multitime correlation functions. There are no multiple points with prob- 
ability one for t > 0. Then we assume that there is a function ]K(s, x; t, y), which is 
continuous with respect to {x,y) G for any fixed G [0, oo)^, such that 



p(t^,x^^^-...-tMiX^^'A= det 



l<j<Nm,l<k<N„ 
l<m,n<M 



for any integer M > 1, any sequence {Nm)m=i of positive integers, and any time 
sequence Q < ti < ■ ■ ■ < Im < oo. Let T = {ti, . . . .tu}- We note that S'^ = 
'YliteT ^t®'^{i) is a determinantal (Fermion) point process on T x M with an operator 
/C given by 



/C/(s,x) = J] I dyK{s,x-t,y)f{t,y), f{t,-)eCo{R),t 



G T. 



When /C is symmetric, Soshnikov [21] and Shirai and Takahashi [23] gave sufficient 
conditions for IK to be a correlation kernel of a determinantal point process. Though 
such conditions are not known for asymmetric variety of processes, which 

are determinantal with asymmetric correlation kernels, have been studied. See, for 
example, [221 112] • If there exists a function K, which has the above properties and 
determines the finite dimensional distributions of the process (S(t),P^), we say the 
process (S(t),P^) is determinantal with the correlation kernelK. \13\ . 

For N E N, the determinant of an x matrix M = {mjk)i<j,k<N is defined 

N 

by sgn(cr) JJ^ mjo-(j), where sgn(cr) denotes the sign of permutation a. Any 

permutation a consists of exclusive cycles. If we write each cyclic permutation as 

a b ■ ■ ■ uj 
b c ■ ■ ■ a 

and the number of cyclic permutations in a given a as ^(cr), then the determinant 
of M is expressed as 

det M = ^ sgn(a) JJ (mabrrihc ■ ■ ■ m^a^ . 
o-e^iv Cj:i<j<e(cT) 

It implies that, with given ai, 02, . . . , oat, even if each element rrijk of the matrix M 
is replaced by rrijk x (aj/af^), the value of determinant is not changed. The above 
observation will lead to the following lemma. 
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Lemma 2.1 Let (H(t),P) and (S(t),P) be the processes, which are determinantal 
with correlation kernels IK and K, respectively. If there is a function G{s, x), which 
is continuous with respect to x eM for any fixed s G [0, oo), such that 

ns,x;t,y) = ^^K{s,x;t,y), s,tG[0,oo), x,y eR, (2.2) 

then 

(S(t),P) = (S(t),P) (2.3) 
in the sense of finite dimensional distributions. 

In literatures, (12.21) is called the gauge transformation and (12.31) is said to be the 
gauge invariance of the determinantal processes. 

2.2 The Weierstrass canonical product and entire functions 

For e Tlo, ^^(M) = iV < oo, with p G No = N U {0} we consider the product 



where 



1 — u if p = 



G{u,p) 



[1 — u) exp 



2 pi (2.4) 



u + — + ■■■ + — 

2 p 



if p G N. 



The functions G{u,p) are called the Weierstrass primary factors. With a > we 
put 



l/a 



{0}- 

For ^ eTlo with ^(M) = oo, we write M«(^, L) for M„(^ n [-L, L]),L> 0, and put 
Ma{C,) = lim Ma(^, L), if the limit finitely exists. If Mp+i(^) < oo for some p G No, 



L— >oo 

the limit 



Up{^,z)= lim Up{^n[-L,L],z)= n g(-,p), zeC (2.5) 



a;G^n{0}<^ 



finitely exists. This infinite product is called the Weierstrass canonical product of 
genus p [L6\. The Hadamard theorem (TB] claims that any entire function / of finite 
order pf < oo can be represented by 

/(z)=;^™e^'(^)n,(e/,z), (2.6) 
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where p is a nonnegative integer less than or equal to pj, Pq{z) is a polynomial 
in z of degree q < Pf, m is the multiplicity of the root at the origin, and ^/ = 
Exe/-i(o)n{o}<= ^x- We give two examples ; 



sin(7r2;) = tczIIo{^z, z), 
Ai(^) = e^'^^'^^U.iU, z) 



(2.7) 
(2.8) 



with ffTTTD . (dH) and 

do -- 

di - 



logAi(O) = - log (^32/3r(2/3) 
Ai'(O) 3V3r(2/3) 35/6(r(2/3))2 



Ai(0) r(l/3) 
For G Tto with ^^(M) = we put 



2tt 



(2.9) 



<l>p(e^,a,;2)=np(r_,e'^,^-a)= J] G 



2; — a 



X — a 



p , a,;zGC. (2.10) 



With ffTT^ we set 



no(e 



u,ie\z), zee, (2.11) 



^Aie\a,z) = $^(r_„e'\z-a 



e'ii{--«)exp 



2; — a 



X 



$o(r,a,z), a,zGC. 



(2.12) 



Lemma 2.2 Lei G with ^^( 
a G supp z ^ a, 



where $^(-,a) = 9<l>^(-, 2;)/92; 



N < 00 and ^^({0}) = 0. Then for 



(2.13) 



$^(e^,a,^) 





1 <^>^(e 




z-a<i.:,(e 




z=a 






a) = {x - 


- 2;)/(a; — a) - 


= {l-z/x)/{l 


e'^'' exp 


[ ^^Aid^) 


lioie -6a,z) 




.Jr. 


no(^^-5a,a) 



(2.14) 
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where the numerator is equal to $^(^^, z)/(l — z/a). From (12. lip , we have 
d 



dz 



+e 



I -^^(rfx)exp 



-cm 



+e'^^^ exp 



Since a G supp is assumed, 

,a) = e^^^exp 



X 



-6a, a). 



It imphes that the denominator of fl2.14p is equal to — a$^(^^,a). Then fl2.13p is 
obtained. ■ 



For G mto with ^^(R) = iV we put 

MA{e)=l -{eA-e){dx). 



(2.15) 



For ^ G mto with ^(M) = oo we write M^(^,L) for Mj^{i n > 0, 

and put M4(^) = lim M_4(^,L), if the limit finitely exists. For ^ G OJIq, p G 

L— >oo 

No, a G R, and 2 G C we define $p(^, a, 2;) = lim $p(^ fl [a — L, a + L], a, 2;) and 

L— >oo 

$^(^, a, z) = lim $^(.^ n [a — L, a + L], a, z), if the limits finitely exist. We note 

L^oo 

that $p(^,a, z) finitely exists and is not identically 0, if Mp+i(r_a,^) < 00, and 
$^(^,a, z) does and $^(,^,0,2) ^ 0, if |Af4(r_aOI < 00 and M2(r_„^) < 00. For 
^ G 9Jto, a G supp ^, the following equalities will hold, if all the entries of them 
finitely exist ; 

z\?({0}) a 



$o(e, «, ^) = no(e, ^)$o(^ n {or, a, 0) (^) 

$o(e n {or, a, 0) = no(e n {-ar , -a)^o{i^'^ n {0}^ a^ 0)2i-««-"», 



a — z 



and then 



$i(e,a,^) = e^(«''^'^)ni(e,^)ni(^n{-ar,- 



x^.„«<=>n{o}^a^o)(i)«"°"^, 

Va/ a — z 



(2.16) 



where 



S{C,a,z)= [ ^^^dx)- [ -adx)+ f -adx). (2.17) 

J{aY X — a J^Qjc X J{0-aY X 
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Lemma 2.3 For a E A, z a 

1 Ai{z) 



^AiUa,z). (2.18) 



z — a Ai'(a) 
Proof. By dHHD and the definition (ETTj) . 

Ai(z) = e'^«$^(a^), ^eC. (2.19) 
As approximations of the Airy function we introduce functions 

Ai;v(2) = e'^o+'^^"JJ( l--)e"/'^^ NeN, (2.20) 
i=i ^ 

where > ai > ■ ■ ■ > are the first zeros of Ai(2;). Since = J2f=i satisfies 

the condition of Lemma [221 (EH with (l2A9ll and Ai'^(aj) = e'^«<l>^(^_^, a^) gives 

' -$^(el,a„^), l<J<Ar. 



2^ flj 



Taking A^ oo, we have (12181) . | 
2.3 Statement of results 

For the solution X{t) = (Xi(t), X2(t), . . . , X7v(t)) of Dyson's model dH]) with /3 = 2 
with the initial state ^(0) = x, we denote the distribution of the process S(t) = 
SjLi '^x^(t) by with = In |T3] we proved that Dyson's model (II. ip 

with (3 = 2 starting from any fixed configuration G OJl is determinantal with the 
correlation kernel given by 

{s,x;t,yj = -= (p az 



27rv-l /r(c^) V -1 

xpsi^{s,x-z) — - — TT (l-—. — -] Psin{-t,w -y) 
w — z -^-^ \ X — z J 

-l{s>t)p,i^{s-t,y-x), (2.21) 

where r(^^) is a closed contour on the complex plane C encircling the points in 
supp on M once in the positive direction, and psin is given by (11.201) . If G TIq, 
by performing the Cauchy integrals (I2.2ip is written as 

K^"" {s,x;t,y) = [ ^^{dx') [ -^^p^^^{s,x - x)^o{i^ ,x' ,y')p^ni{,-t,y' - y) 



-l{s>t)p,i^{s-t,x-y). (2.22) 



Then the following is obtained for the process (S^(t), ) with H_4(t) = X]j=i ^Yj{t), 
where Y{t) = (Fi(t), . . . , YNit)) is given by (|rT2il . 
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[s, X- t,y) = I ^''{dx') I q{0, s,x- x' , y')q{t, 0, y' - y) 



Proposition 2.4 The process (S_4(i(:),P^ ), starting from any fixed configuration 

G OJlo with ^^(M) = N < oo, is determinantal with the correlation kernel 
given by 

dl_ 

-IR V— 1 

-lis>t)qit,s,x-y), (2.23) 
where q is given by lil.30\) . 

We introduce the following conditions: 

(C.l) there exists Co > such that |Af4(^)| < Co, 

(C.2) (i) there exist a e (3/2,2) and Ci > such that M^iO < Ci, 
(ii) there exist /3 > and C2 > such that 

Miir^a^^^"^^) <C2{\a\yi)-^ for all a G supp 

We denote by X-^ the set of configurations ^ satisfying the conditions (C.l) and 
(C.2), and put = X-^ n Mq. 

Theorem 2.5 IfC, G Xq-, the process (S_4(t),P^) is well defined with the correlation 
kernel 

f f dy' 

^A{s,x;t,y) = / ^{dx') / -j= q{0, s,x - x')<i>A{^,x' ,y')q{t,0,y' - y) 

-lis>t)qit,s,x-y). (2.24) 

In the proof of this theorem, a useful estimate of $_4 in fl2.24p is obtained (Lemma 
14.31 (ii)). By virtue of it, we can see 

K^^{t,x;t,y)K^^(t,y;t,x)dxdy ^ ^{dx)l{x = y) as t — > (2.25) 

in the vague topology. Then Theorem 12.51 gives an infinite particle system starting 
form the configuration ^. 

The main result of the present paper is the following. 

Theorem 2.6 (i) Let (■) be the configuration Iil.l5\) . Then 

(S^(t),P^^)^(S^(t),PAi) as iV-.oo 

in the sense of finite dimensional distributions. Here the process (H_4(t),PAi) is de- 
terminantal with the correlation kernel U.23\) with setting U.26\) and pf{s,x]t,y) = 
PAi(t — s,y\x) with U.31\) . that is 

sr^ f dz 1 Ai(z) 

KAi{s,x;t,y) = > / ^=pAi{s,x\a) ——pAi{-t,z\y) 

T-^,^Jy/^RV^ z-aAi{a) 



aeAi-i(O) 

-l{s>t)pAi{s-t,x\y). (2.26) 
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(ii) Let (S_4(t),PAi) be the process, which is determinantal with the extended Airy 
kernel 17371 ). Then 

(S^(t + ^^),PAi)^(H^(t),PAi) as (2.27) 
weakly in the sense of finite dimensional distributions. 

2.4 Remarks on Extensions of the Results 

(1) By definition (12.151) . A'/4(^^) = 0. The asymptotic property of the zeros (11. lip 
imphes 

C>)^ (M,(a))" = $^^<oo, if a>p/ = | (2.28) 

In general, order of growth pj of a canonical produce (12. 5p is equal to the convergence 
exponent pj of the sequence of its zeros [16j. For Ai(2;), pj = 3/2. The function 
C'^{a) may be called the Airy zeta function [30j, which is meromorphic in the whole 
of C [7j. Moreover, we know 

C-^(2) = Mi(ei?>) = E i = ^? < ^ (2.29) 

aeA 

with (12. 9p . Then satisfies the conditions (C.l) and (C.2) : G X"^. Since 
^ 2Jto, Theorem 12. 51 guarantees the well-definedness of the infinite particle system 
(S_4(t),P^). (Its equivalence with (S(t),PAi) is stated in Theorem 12.61 (i).) Note 
that the negative divergence (11.161) of the drift term Djs^j^t of (I1.12p m N oo for 
t < oo corresponds to that C'^(l) = — Sagyi(l/ct) = oo. This fact and (12.291) mean 
that the Airy function has genus 1 [HI 130] . 

Examples of infinite particle configurations in other than ,^_4 are given as 
follows. For K > 0, we put 

g'^ix) = sgn(x)|x|^ xER, and ri%-) = J^^gHni')- 

eez 

For any k > 1/2 we can confirm by simple calculation that any configuration C, G 'OJIq 
with supp ^ C supp rj'^ = {g'^{i) : i E Z} satisfies (C.2)(i) with any a E {1/k,2) 
and some Ci = Ci{a) > depending on a and (C.2)(ii) with any j3 E (0,2k — 1) 
and some C2 = C2{P) > depending on /3. Assume that ^ E VJIq is chosen so that 
supp ^ C supp rj'^ for some k > 1/2 and |M4(^)| < 00. Then ^ E X^. The fact 
(II. lip implies that this assumption can be satisfied only if k G (1/2,2/3]. 

(2) If there exists, however, f3' < {/3 - 1) A {13/2) for i E such that "^{x E i : 
i{\x — \x\^ , X + \x\^ ]) > 2} = cxD, then ^ does not satisfy the condition (C.2) (ii). In 
order to include such initial configurations as well as those with multiple points in 
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our study of the process (E!_4(t), P^) with = oo,we introduce another condition 
for configurations: 

(C.3) there exists n G (1/2,2/3] and m G N such that 



We denote by the set of configurations ^ satisfying (C.l) and (C.3) with 

K G (1/2, 2/3] and m G N, and put 



2)-" = U U ?^ 

Ke(l/2,2/3] meN 



A 



Noting that the set G 9Jt : m(^, k) < m} is relatively compact for each k, G 
(1/2, 2/3] and m G N, we see that 2)-^ is locally compact. 

In the present paper, we report our study of the relaxation process (S(t),PAi) 
from a special initial configuration to the stationary state fiAi- We expect that ^Ai 
is an attractor in the configuration space 2)"^ and is a point included in the basin. 
Motivated by such consideration, we are interested in the continuity of the process 
with respect to initial configuration. We have found, however, that if ^(M) = oo, the 
weak convergence of processes in the sense of finite dimensional distributions can not 
be concluded from the convergence of initial configurations in the vague topology. 
Following the idea given by our previous paper [13], we introduce a stronger topology 
for 2)-^. 

Suppose that ^,^„ G 2)-^,n G N. We say that ^„ converges ^_^-moderately to ^, 

if 



lim $yi(^n, V — 1, ■) = '^'^(^) V — 1, ■) uniformly on any compact set of C. (2.30) 

n— >oo 

It is easy to see that fl2.30l) is satisfied, if the following two conditions hold: 



lim sup lim 

L— >oo n>0 M—*oo 



M^(e„,M)-M^(e„,L) 



0, 



lim sup 



(2.31) 



(2.32) 



By the similar argument given in [IB], the following statements are proved. 

(i) If ^ G 2)-^, the process (S_4,P^) is well defined. 

(ii) Suppose that e 2);^„,?i G N, for some k G (1/2,2/3] and m G N. If C.n 
converges $^-moderately to ^, then (S_4,P^) (S_4,P^) weakly in the sense of 
finite dimensional distributions as n ^ oo in the vague topology. 

Moreover, we can show /iAi(2)'^) = 1. By this fact and the above mentioned 
continuity with respect to initial configurations, we can prove that the stationary 
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process (S_4(t), Pai), which is determinantal with the extended Airy kernel (11.371) . 
is Markovian [14J. 

(3) As mentioned in Introduction, the purpose of the present paper is to give a 
method for asymmetric initial configurations to construct infinite particle systems of 
Brownian motions interacting through pair force 1/x. In order to clarify the results, 
we have concentrated on the case in this paper such that the initial configuration 
is (Theorem 12. 6p or its modification ^ G Xjf; see the condition (C.l) with 
(I2.15P (Theorem 12.51) . In the former case the constructed infinite particle system 
(S_4()f:), PAi) has the stationary measure /^Ai, which is obtained in the soft-edge scaling 
limit of the eigenvalue distribution in GUE well-studied in the random matrix theory. 
Thus we have specified the entire function used in our analysis in the form $^(^, z) 
given by (12. lip , which is suitable for the Airy function (see (12.190 ). The point of 
our method is to put the relationship between the entire function appearing in the 
correlation kernel K^, the "typical" initial configuration C,f, and the drift term in the 
SDEs providing finite-particle approximations. By the same argument as reported 
here, the following will be proved. Let / be the entire function such that /(O) ^ 0, 
it is expressed by the Weierstrass canonical product of genus one, ni(^/, z), and the 
zeros can be labelled as < |xi| < \x2\ < ■ ■ ■ . Then with ^/ = J2'jLi we put 



^ 1 

= 1 



and introduce the A^-particle system 

Y/{t)=X,{t) + D{^f,N)t, 1<J<N, te[0,oo), 

where X(t) = {Xi(t), . . . , Xiy{t)) is the solution of Dyson's model (II. ip starting from 
the first N zeros of /, Xj{0) = Xj,l < j < N. Then Sj(t) = ^Y^'{t) converges 
to the dynamics in — oo, in the sense of finite dimensional distribution, which is 
determinantal with the correlation kernel 

¥^/{s,x;t,y) = ^f{dx') -j=psi^{s,x - x')^i{Cf,x',y')psin{-t,y' - y) 

Jm. Jy^R y—i 

-l{s > t)psin{s -t,x - y). (2.33) 

Moreover, even if the initial configuration is different from C,f, but it satisfies the 
condition 



[ -iif-Oidx) 

J-L X 



for any L > with a positive finite Cq independent of L, then the process starting 
from ^, is well-defined. In general, the obtained dynamics with an infinite number 
of particles is not stationary, while Theorem 12.61 gave the example which converges 
to a stationary dynamics (S_4(t),PAi) in the long-term limit. 
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3 Properties of the Airy Functions 
3.1 Integrals 

By the fact Ai"(x) = xAi(x), the following primitive is obtained for c 7^ [30\ . 

du (Ai(c(u + x))f = {u + x)(Ai(c(M + x))f - ^(Ai'(c(u + (3.1) 



du Ai(c('U + x))Ai(c(-u + y)) 

Ai'(c(M + a;))Ai(c(u + y)) - Ki{c{u + x))ki{c{u + y)) 



c^{x-y) 

By setting c = 1 and integral interval be [0, 00) in (13.21) . we obtain the integral 

M{x)M'{y) - M{x)M{y) 



(3.2) 



du M{u + x)PA{u + y) 



x-y 



since lim^^-^oo Ai(a;) = lim^^^oo Ai'(x) = by (II. 5p . If we set ?/ = a G ^ and 
X = z ^ a, then 

r , ^-f ^ Ai{z)Ai'{a) 

/ du Ai(u + z)Ai(u + a) = — — 

Jo z-a 

since Ai(a) = 0. Then we have the expression 



1 Ai{z) 



POO 

/ du Ai{u + z)Ai{u + a) (3.3) 
Jo 



z-aAi'{a) (Ai'(a))2 J^ 
for a E A, z a. 

3.2 Airy transform 

The following integral formulas are proved [11] . 

Lemma 3.1 For c > 0,x,y 

[ du e'="Ai(M + x)Aiiu + y) = ^Le-(-s/)V(4c)-c(x+y)/2+c3/i2^ ^3 4^ 

Jr V 4:710 

dy [ du e™Ai(u + x)Ai{u + y) = e'^^+'^'/l ^3 5) 

: Jr 

Proof. Consider the integral 

/ = / due'^'^Ailu + x)Ai{u + y) 
Jr 

POO pO 

= du e™Ai(M + x)Ai{u + y) + du e''''Ai{u + x)Ai{u + y). 
Jo J -00 
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By the definition of the Airy function (11. 4p . for any 1] > 0, 
1 



+ 



du 



du 



+ 



— oo 

1 

1 

1 

1 

(2^ 
1 



(27r)2 

(27r)2 
dz 



dz 



dw gV^i-^ /3+{"+^02W/3+(2/+u)m>} 



32=r; J<sw>c—rj 

dz 



^^^V^{z^ /3+{-^+x)z+'W^ /3+{u+y)w} 
i3z=T] J'^wKc—rj 

3 /Q^^^_i_.,,,3 / 



dw e^-T{23/3+a;2+w3/3+j/w) / ^{c+y^{z+w)}u 



isz=ri J<sw>c—ri 

dz 

'Csz=r] J'Csw<c—-q 

dz 



oo 







1 



dw eV^(-'V3+2:^ W/3+yt«) 

1 



isz=ri J<sw<c—ri 
3 , 



C + V^(-2 + "U^) 



' isz=r) 

It is equal to the integral 
1 



dw — 



<sw<c—ri 



dw 



\sw>c—'ri 



-l{w^/3+yw) 



C + \f—\{z + w) 



2H 



<SZ=T) 



dw ■ 



27rv — 1 Jc w 



-Ic - z) 



where C is a closed contour on C encircling a pole at w = y— Ic — z once in the 
positive direction. By performing the Cauchy integral, we have 



1 

2^ 



dz eV^(^V3+x^)gV^{v^c-^)V3-2/(c+v^z) 



<sz=r] 



]_ c^/3-cy 

2tt 



<sz=r] 



By performing a Gaussian integral, we have fl3.4p . Since 



y — {x — c^) > — cx + 



the Gaussian integration of (13. 4p with respect to y gives (13.51) . i 

By setting c = t/2 > in (13. 5p and c = (t — s)/2 > in (13. 4p . respectively, we 
obtain the equalities 

dy due^^M{u + x)ki{u + y) = g{t,x), 
Jr Jr 

I du e"(*-^)/'Ai(M + a;)Ai(M + y) = ^^^q(s, t,y-x) 
Jr 9{s,x) 
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with fOHD and fOOj) . Thus we have defined g{t,x) for any t G M by ffT:28D and 
PAi{s, X] t, y) for any s, t G M, s 7^ t by fll.311) . As special cases, we have 

PAi{t,y\x) = g{t,y)q{0,t,y - x), (3.6) 

PAii-t,y\x) = — — -q{t,0,y-x), t>0,x,y,eR. (3.7) 
g{t,x} 



If we take the c — hmit in (13. 5p . we obtain 

[ [ dxAi{^-x)Ai{C -x) = 1. 



The expression (13. 4p and the above result implies the orthonormality of the Airy 
function in the sense; 

(^J duAi{u + x)Ai{u + y)^ dy = 5^{dy). (3.8) 
The Airy transform f{x) t— > ^p{^) is then defined by 

V{0= f dxf{x)Ai{^ + x), (3.9) 



and the inverse transform is given by /(x) = d^ (f{^)Ai(^ + x). Now a parameter 

Jr 

c G C is introduced and the family of functions are defined as {wc{x) = Ai{x / c) / \c\} . 
The Airy transform (13. 9p is then generalized as 

M0= [ dUix)w,i^ + x) = ^ [ rfx/(x)Ai^^ + '' 



Lemma 3.2 The Airy transform with c of the normalized Gaussian function f{x) = 
e-^7v^ «s given by cp^iO = \c\~^e^^+^/^^'^''^^/^'^'''^Ai{^/c+l/{lQc'^)). That is, 



Proof. By the definition of the Airy function (II. 4p . 
I = [ dx ^e'"'^ AS. ' ^ ^ ^ 



IT \ C 



— [ rffcev^'^'/^+v^^fc/c 1 f dx exp ( -x^ - V^-x 
2vr y/TT Jr V c 
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By performing the Gaussian integral we have 



1 /• / — A;3 ^ — 



By completing a cube, we find the equality 



3 c 4c2 3 V 4c2 



c 16c4 

By using the definition of the Airy function (11.41) . f l3.10p is obtained, i 
For t > 0,y,u E M. we will obtain the equality 

/ -4^Ai(u + z)q(t,0,z-y)= [ dx ^e'""^ Ai ( \^^2i x + u + y - 

by changing the integral variable as z ^-^ x = {z — y + t'^/A)/y/—2t. If we set 
^ = {u + y- tV4)/y^ and c = l/V^, the RHS is identified with the LHS of 
dXTU]) . Since 



^ 1 1/1 



' 2 24 



ut 



( 13:T011 of Lemma O with ([LSHD gives 



Ai(u + 0, 2; - = g{t, y)e~'''/^Ai{u + y), t>0,y,ue 



Combination with fl3.7p gives 

dz 



-1 



Ai(M + 2)pAi(-t,^|l/) = e~"*/^Ai(M + ?/), t>0,y,ueR (3.11) 



3.3 Fourier-Airy series 

Let us consider the integral 



Iff 



POO 

/ dx Ai{x + ae)Ai{x + aei), ai^a^iEA. 
Jo 



In the case i ^ i', the formula (13. 2p gives 



Ai'(a^)Ai(a^O - Ai(a£)Ai'(a^O 
— — U, 

0,(1 — 0.^' 
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whereas if i = i', the formula (13. ip gives = (Ai'(a^))^. Therefore the functions 

Ai(a; + at 



Ai («£ 



(3.12) 



form an orthogonal basis for / G L^(0, cxd) (see Sect. 4.12 in [27j). The completeness 
of fl3.12p is also established : 



^ Ai(x + af)Ai(y + a^) 



(Ai'(a,))^ 



(3.13) 



Then for any / G L'^{0, oo), we can write the expression 



Ai (a^) 



and call it the Fourier-Airy series expansion. The coefficients q of this expansion 
are determined by q = {j^ dx f{x)Ai{x + a^)}/Ai'(a^). 

4 Proof of Results 



4.1 Proof of Proposition 12.4 

With (fTT^ we put 



g{s,x) = exp <^ ~D_A 



By the definition fll.30p of q, for s > 0,x,x' E M, we have 
Psm(s, (x - Daj^s - s^/A) - x') 



g(0, s,x — x') 
exp 



X X 

4 



exp 



h X + X 



^(s,x)e-^^iv-', 



and for t > 0, G M, we have 

(l{t,0,y' -y) 



9{t, y) 
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Similarly, we have 

Psinjs - t, (x - Dj^^s - gV4) - (y - Dj,^t - tV4)) 
q{t,s,y- x) 

exp -^^A[--y-DAAs-t)-^ 



X 



exp 



-Da, 



[s-t) + 



' 4 



(x-y) 



9{t, y) ' 



Then we have 



K'^ \s,x-DAj^s-—;t,y-DAj-—\ 



9{t, y) 



e{dx') 



dy' 



-II 



xg(0, s,x- x')e"^^iv-'$o(^^, x', y')e''^Ny' q{t^ 0, y' - y) 

-l(s > t)q{t,s,x-y) 



The identity (12.121) implies 



By the gauge invariance of determinantal processes (Lemma 12. II) . the proof is com- 
pleted. I 



4.2 Proof of Theorem [23 



First we prepare some lemmas for proving Theorem 12.51 

Lemma 4.1 Let a G (1,2) and 6 > a — 1. Suppose that Mad) < oo and put 
Lo = Lo(a,<5,0 = {2M4^))^/^'-"+^l Then 

M^{i,L)<L\ L>Lo. 

Since this lemma was proved as Lemma 4.3 in [13j, here we omit the proof. 

Lemma 4.2 If satisfies (C.2) (i) and (ii), for any 9 G (a V (2 — /?), 2) there exists 
C = C{Ci,C2,e) > such that 



f {--—) 

i{o,a}c \x x-aj 



i{dx) 



<C|aVl|^-\ a G supp (e-5o)- (4.i; 
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Proof. We divide {0, aY into three sets Ai = {x E {0, aY : |x| < |a|/2}, A2 = {x E 
{0,aY : \a\/2 < \x\ < 2\a\}, and = {x E {0, a}^ : 2\a\ < \x\}, and put 



1 



tXy ijC Gj 



i{dx) J = 1,2,3. 



When X E Ai, — a,^| > 3a^/4 and |x + a| < 3|a|/2, and then 

lal |x + a\ 



-^dx) < 2Mi U 



By Lemma [4.11 for any 6 > a — 1, we can take C > such that 

h < C|a V 1|^ 
When X E A2, \x + a\ < 3\a\ and |a|/|x| < 2, and then 

lal |x + a' 



A2 



From the condition (C.2) (ii) 

h < eCala V 1|^-^. 
When X E A^, \x — a\ > \x\/2, and then 





\a 




\x\ 


X ■ 


- a\ 



^{dx) < 2"->|°-W„(0 



From the condition (C.2) (i) 



I3 < 2"~^Ci|a V 1 



la-l 



Combining the estimates ( lO) . ( Oil and ( lOll . we have (lO) . 



(4.2) 



(4.3) 



(4.4) 



Lemma 4.3 (i) //^ satisfies the conditions (C.2) (i) and (ii), for any 9 E [aM 
(2-/3), 2) there exists C = C{Ci, C2, 6) > such that 

|$i(e, a, V^y)\ < exp [C{{\yf v 1) + (ja^ V 1)}] , (4.5) 

for y eM and a E supp ^. 

(ii) //^ satisfies the conditions (C.l), (C.2) (i) and (ii) , for any 6 E (a V (2-/3), 2) 
there exists C = C(Co, Ci, C2, 6*) > such that 

|$^(e, a, V^y)\ < exp [C{(|y|^ V 1) + (|a|% 1)}] , (4.6) 

for y eM. and a E supp ^. 
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Proof. We prove (i) of this lemma. From the condition (C.l) and the relation 
fl2.12p . (ii) is easily derived from (i). We first consider the case that a = G supp ^. 
Remind that 



$1(^,0,2;) = Hi (e, 2;) =exp 
When 2\z\ < \x\, by using the expansion 



J{o}c I- V xJ x) 



log 1 



X 



keN 



1 



k \x 



we have 



Then 



log 




z 




z 




+ - 


< 








X 




X 



|ni(en[-2|z|,2|^|]^;2)| <exp(|^|2 \ <exp 

L J\x\>2\z\ X ) 



'\x\>2\z\ 

On the other hand |1 — < e'^'/'^L Then 



i;2rM„(er . (4.7) 



|ni(en[-2|^|,2|z|],z)| <exp(2|z| [ ^f{dx)]=exp{2\z\M^{^,2\z\)]. 

L J\x\<2\z\ \X\ ) ^ 

(4.8) 

From (14.71) and (14.81) . with the condition (C.2) (i) and Lemma [4.11 we see that for 
any 6* G (a V (2 - /?), 2) there exists C = C{Ci, 9) > such that 



|$i(e,0,z)| <exp [C{i\zfyi)}] 



(4.9) 



for z G C and a G supp ^. 

Next we consider the case that a G supp C, and a 7^ 0. By the conditions (C.l) 
and (C.2) the equality (EZTE]) is valid. By 

|ni(e, z)U,{^ n {-ar, -a)\ < exp [c{{\zf V 1) + {\af V 1)}" . 

By the condition (C.2) (ii) 

I$o(e^'^n{0r,a2,0)| <exp{|apMi(r_,.e^2>)| <^^p|^^(|^|^^^2-/3|^ 

and \{y/—ly/a)^^^^^^a/ (a — y/—ly)\ < 1. Now we evaluate S{^,a,z). 

z — a z — a 



S{i,a,z) 



{z — a) 



^{dx) 



+^ ^^({0})-- + --— m-a}) 

a a a —a 



,z — a) 



I (^-i)ew-^ii« + i + «{o}) + «{-a}). 

J{Q,aY \x-a xJ a 
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From Lemma [4.21 and the fact < C2 and then \2z/a\ < 2a/C2^|z|, we have 
\S{^, a,z)\< C\z - a||a V ll^^^ + 2v^|;z| + 3 < C"|(|y|^ V 1) + {\af V 1) 
for some C > 0. This completes the proof. 1 

Proof of TheoremlEE Note that ^ n [-L, L], L > and ^ satisfy (C.l) and (C.2) 
with the same constants Co, Ci, C2 and indices a, (3. By virtue of Lemma 1^13) (ii) we 
see that there exists C > such that 



\<^M n [-L, L], a, v^y)| < exp \c[{\y\ V 1)^ + (|a| V 1)^ 
VL > 0, Va G supp Vy G M. Since for any y G M 

we can apply Lebesgue 's convergence theorem to fl2.23p and obtain 

lim K^"'"'^'-^^ (s, x; t, y) = (s, x; t, y) . 



L— >oo 



Since for any (s, t) G (0, 00)^ and any finite interval I C 



sup 

x,y(il 



«n[-L,L] 
.4 



's,x;t,y) 



< 00, 



we can obtain the convergence of generating functions for multitime correlation 
functions g^^^-^'^^[x] ^ G^x] as L ^ 00. It implies ^ Pj^ as L cx) 

in the sense of finite dimensional distributions. Then the proof is completed. | 



4.3 Proof of Theorem [276 



(i) It is clear that C,a is an element of j£^. (See the item (1) of Sect. 2.4.) Then 

by Theorem [23] (2^ (t),P^) ^ (H^(t),P$f) as iV ^ 00 in the sense of finite di- 
mensional distributions, where (S_4(t), P^) is the determinantal with the correlation 
kernel 

IK^(s,a;;t,?/) = I ^^(rfa) /" g(0, s, x - a)$^(^^, a, 2;)g(t, 0, 2; - y) 

> t)q{t,s,x - y). 

Using the equalities (13.61) and (13. 7p and the definition (11.311) . we have 
jirUt + \ f ^ /J \ f PA\{s,x\a) 1 M{z 



-1 (7(5,3;) 2; — aAi'(a^ 



9{'t^y)PAi{-t^z\y) 



-l(^>^)4^PAi(5-t,x|i/) 



^(^ 
^(s,x) 



5((s,x)' 
]KAi(s,x;t,2/) 
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with (12.261) . where we have used fl2.18p of Lemma 12.31 By the gauge invariance, 
Lemma [^TT| (S_4(t),P^) = (S_4(t),PAi) in the sense of finite dimensional distribu- 
tions. 

(ii) If we use the expression (13. 3p . (12.261) becomes 



KAii-s, x;t,y) = / -j=PAiis, x\a) 

. .-1 .„^ Jx/^R V — 1 



aeAi-i(O) 
X 



1 

— — / duAi{u + z)Ai{u + a)pAi{-t,z\y) 
1 [aW Jo 



(A: 

> t)pAiis - t,x\y). 
By (13. lip , the first term of the RHS equals 

1 

PAi{s,x\a) / rfue""*/2Ai(u + y)Ai(M + a). 

Since s > 0, we can use the expression (I1.27P for PAi{s, x\a) and the above is written 
as 

Hdu [ dw e--/^+-/^Ai(. + y)Aiiw + x) V ^'^""^^^ft. ^ "'^ • 
Jo Jr (Ai(a^))^ 

From the completeness (I3.13p . the above gives 

]KAi(s, x; t, y) = KAi(t - s, y\x) + R{s, x; t, y) 
with the extended Airy kernel Kai given by (I1.37P and 



R{s,x;t,y) = j du dw e-'''/^+'"'/^Ai{u + y)Ai{w + x) 

J —oo 

Ai{u + ai)Ai{w + a^) 



X 



(Ki(a,))' 



Since for any fixed s, t > limg^oo \R{s + 6,x]t + 6,y)\ — > uniformly on any 
compact subset of M^, (I1.36P holds in the same sense. Hence we obtain (I2.27p . This 
completes the proof, i 
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